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^D I We present a direct measurement of the parity violation parameter Af,, de- 

^ ' rived from the left-right forward-backward asymmetry of b quarks tagged via 

leptons from semileptonic decays. The lepton identification algorithm com- 



Qh! bines information from tracking, calorimetry and from the SLD Cherenkov 

Ring Imaging Detector. The value of Ai, is extracted using a maximum like- 



lihood fit to the differential cross section for fermion production. Vertexing 
information and decay kinematics have been used to discriminate among the 
kS ' different sources of tagged leptons. A new treatment of mixing effects and 

C^ ' of background contamination has been introduced and a new vertexing al- 

gorithm has been used in the muon analysis. Based on the 1993-1998 SLD 
sample of 550K hadronic Z^ decays with highly polarized electron beams, we 
have measured Af, with a ~ 3% statistical error. 
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1 Introduction 



Parity violation in the Zff coupling can be measured via the observables Af = 
2vfaf/{v'j + Of), where Vf and a/ represent the vector and axial vector couplings to 
fermion /. The Born-level differential cross section for the reaction e~^e^ ^ Z^ ^ f f 
is 

dcTf / rfz OC (1 - AePe){l + Z^) + 2Af{A, - P,)z , (1) 

where Pg is the longitudinal polarization of the electron beam (Pg > for right- 
handed (R) polarization) and z = cos 6 is the direction of the outgoing fermion 
relative to the incident electron. In the presence of e~ beam polarization, it is 
possible to construct the left-right forward-backward asymmetry 



^/ (. _ [4{^) - 4{-^)] - [4(^) - 4(-^)] _ \Pe\Af2z 
'^"^ ' [4{^) + crii-z)] + [aiiz) + aU-z)] (1 + ^^) ' ^ ^ 
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for which the dependence on the initial state coupling parameter Af> disappears, 
allowing a direct measurement of the final state coupling parameters Af. Thus elec- 
tron beam polarization permits a unique measurement of Af, independent of that 
inferred from the unpolarized forward-backward asymmetry |I[] which measures the 
combination Af.Af. In addition, the quantity Afy is largely independent of propaga- 
tor effects that modify the effective weak mixing angle, and so is complementary to 
other electroweak measurements performed at the Z^ pole. In particular the Stan- 
dard Model expectation Af, = 0.935 has only a very slight dependence on the top 
quark and Higgs boson masses. 

In this paper we present a direct measurement of Ab based on identified leptons 
from semileptonic B hadron decays. The analysis is based on the full 1993-1998 
SLD data sample of 550,000 Z^ decays and presents the improvements obtained 
with the addition of vertexing information provided by the new vertex detector 
(VXD3) installed in 1996. The measurement complements other direct measure- 
ments of Af, performed at SLD, that use momentum- weighted track charge [Q, ver- 
tex charge and identified kaons @] to determine the sign of the underlying quark 
in bb events. The lepton total and transverse momentum (with respect to the near- 
est jet), the mass of the event and some topological decay information are used to 
classify each event by deriving probabilities for the decays (Z° — > 66, 6 — > lepton), 
(Z° ^ bb,b —>■ c ^ lepton), (Z° ^ bb,b -^ c ^ lepton), (Z° —* cc,c ^ lepton), and 
(Z° -^ background) f\ The lepton charge (Q) provides quark-antiquark discrim- 
ination, while the jet nearest in direction to the lepton approximates the quark 
direction. The parameter A^ is then extracted by a maximum likelihood fit of these 
data to the polarized differential cross section, taking into account the effects of 
hard gluon radiation. Although in this approach the polarized asymmetry (2) is 
not explicitly formed, the result for Ab maintain its insensitivity to the initial state 
couplings. 



"'^ leptons from light hadron decays, photon conversions and misidentified leptons 



2 Data Selection and Lepton Identification 



The SLAC Linear Collider and its operation with a polarized electron beam have 
been described in detail elsewhere [§. During the running period from 1993-98, the 
SLC Large Detector (SLD) recorded an integrated luminosity of 19.1 pb~^ with a 
luminosity-weighted electron beam polarization of \Pe\ = 0.729 ±0.004 (1997-98) at 
a mean center of mass energy of 91.27 GeV. 



Charged particle tracks are reconstructed in the Central Drift Chamber |T0 
and the CCD-based vertex detector [|rT|], in a uniform axial magnetic field of 0.6T. 
The combined momentum resolution in the plane perpendicular to the beam axis is 



5p±/p± = V(-Ol)^ + (-0026 p^/GeVy . 

The Liquid Argon Calorimeter (LAC) [|12| measures the energies of charged and 
neutral particles and is also used for electron identification. The LAC is segmented 
into projective towers with separate electromagnetic and hadronic sections. In the 
barrel LAC, which covers the angular range |cos^| < 0.82, the electromagnetic 
towers have transverse size ~ (36 mrad)^ and are divided longitudinally into a front 
section of 6 radiation lengths and a back section of 15 radiation lengths. The barrel 
LAC electromagnetic energy resolution is as/ E = 1^%/ J E{GeV). 

Muon tracking is provided by the Warm Iron Calorimeter (WIC) [|l3l . The WIC 
is 4 interaction lengths thick and surrounds the 2.8 + 0.7 interaction lengths of the 
LAC and SLD magnet coil. Sixteen layers of plastic streamer tubes interleaved with 
2 inch thick plates of iron absorber provide muon hit resolutions of 0.4 cm and 2.0 
cm in the azimuthal and axial directions respectively. 



The Cerenkov Ring Imaging Detector (CRID) |jT^ measures the velocities of 
charged tracks using the angles of Cerenkov photons emitted in liquid and gaseous 
radiators. Only the gas information has been included in this analysis, since the 
liquid covers only marginally the interesting momentum region (p > 2 GeV/c). 
Electrons are distinguishable from pions in the region between 2 and 5 GeV and 
the muon identification (because of pion rejection) also improves considerably in 
this region. Kaon and proton rejection also helps the muon identification up to 
momenta of 15 GeV. 

Hadronic events are selected by requiring at least 15 GeV of energy in the LAC 
and at least six tracks with p_|_ > 250 MeV. Approximately 550,000 events were 
found in the 1993-98 data sample, with negligible background. Jets are formed 
by combining calorimeter energy clusters according to the JADE algorithm |l^ 



with parameter ycut = 0.005. The jet axis closely approximates the 6-quark di- 
rection in Z*^ —^ bb events, with an angular resolution of ~ 30 mrad. An electron 
or muon tag is used to select semi-leptonic decays. Electrons are identified with 
both LAC and CRID information for tracks with p > 2 GeV in the angular range 
I cos^l < 0.72. Calorimeter information is used to build discriminant variables which 



exploit the characteristics of electromagnetic showers, including transverse and lon- 
gitudinal shower development shapes, and matching of LAC energy and track mo- 
mentum. The GRID information is stored in likelihood functions corresponding to 



each particle type hypothesis |T5[. These quantities are used as input variables to 



a Neural Network, trained on Monte Carlo tracks [0. The effficiency (purity) for 
electron identification is on average 62% (70%) and over 78% (80%) for electrons 
with momenta greater than 15 GeV/c. This estimate includes electrons from pho- 
ton conversions as signal. As pion misidentification is the largest contribution to 
the background, the simulation has been verified using charged pions from recon- 
structed K^ —>■ 7r"'"7r~ decays. The fraction of such pions misidentified as electrons 
is (1.23 ± 0.15)%, consistent with a MC expectation of (1.36 ± 0.07)%. Electrons 
from photon conversions are removed from the sample with a 70% efficiency. The 
remaining photon conversion background is clustered at low momentum, away from 
most of the signal region. 

Muon identification is performed for tracks with p > 2 GeV in the angular 
range | cos6\ < 0.70, although the muon identificaton efficiency falls off rapidly for 
I cos^l > 0.60 (in the region between the barrel and the endcaps). CDC tracks are 
extrapolated along with the associated error matrices, including multiple scattering, 
and matched with hit patterns in the WIC. For | cos^| < 0.60, 87% of the simulated 
muon tracks have successful matching between the CDC and the WIC. The second 
step of the muon identification exploits the information from the CRID. The CRID 
k — fi separation variable alone rejects 51% of the remaining k and p (with only 
2% loss in the signal), while, for p < 6 GeV, the tt — /i separation variable rejects 
37% of TT (with 5% loss in the signal). Since the CRID information is intrinsically 
momentum dependent, different sets of cuts on the distributions of the discriminant 
variables have been optimized in different momentum regions to achieve best purity 
and efficiency. The purity of the final sample is improved by requiring that the 
candidate muons fully penetrate the WIC, and by applying further cuts on the 
number of hits associated with the tracks, on the x^ of the CDC/WIG matching 
and on the x^ of the fit of the track in the WIC. MC studies show that pion punch 
through background is negligible. Muons from pion and kaon decays and hadronic 
showers are a significant background, but fall off rapidly with increasing momentum. 
From a study on a pure pions data sample, obtained from kinematically selected 
K^ -^ TT+TT" decays, ~ .3% of pions, with p > 2 GeV, are identified as muons. The 
muon identification efficiency is over 81% with a purity of 68% (8% misidentified 
tracks and 24% muons from light hadron decays) for p > 2 GeV and | cos^| < 0.60. 



3 Monte Carlo Simulation 



The likelihood that an identified lepton comes from each of the physics sources {b -^ 
l,b ^ c —>■ I, c ^ I, background etc.) relies directly on MC simulation. Z^ decays are 
generated by the JETSET 7.4 program ||18i. The B hadron decay model was tuned 



to reproduce existing data from other experiments. Semileptonic decays of B mesons 
are generated according to the ISGW formahsm [jl9| with a 23% D** fraction, while 
semileptonic decays of D mesons are generated with JETSET with the 1994 Particle 
Data Group branching ratios |^0|]. Particularly important experimental constraints 
are provided by the B -^ lepton and B ^ D momentum spectra measured by CLEO 
21[] [ p^ , the D — >■ lepton momentum spectrum measured by DELCO [^, and the 



B -^ hadron multiplicities measured by ARGUS |24 



The SLD detector response is simulated in detail using GEANT ||2^ and has 
been checked extensively against Z^ data. 



4 Vertex Mass Reconstruction 



Vertex identification is done topologically, by searching for space points in 3D where 
track density functions overlap 0. Each track is parametrized by a Gaussian prob- 
ability density tube with a width equal to the uncertainty in the measured track 
position at the IP. Points that are characterized by a large overlap of these Gaussian 
probability tubes are considered candidate [seed) vertices. By clustering maxima in 
the density distribution, secondary vertices are found for the two hemispheres. The 
efficiency for reconstructing a vertex in the same hemisphere as the lepton is ~ 66% 
(1996-98). The mass of the secondary vertex is calculated using the tracks attached 
to the vertex itself. Each track is assigned the mass of a charged pion and the in- 
variant mass of the vertex is thus calculated. This is then corrected to account for 
neutral particles by using kinematic information. By comparing the vertex flight 
path and the momentum sum of the tracks associated to the secondary vertex, one 
calculates a minimum amount of missing transverse momentum to be added to the 
invariant (raw) mass. This is done by assuming that the true quark momentum is 
aligned with the flight direction of the vertex. The so-called P^-corrected mass is 
then given by: 



MvTX = ^/M^^^ + P,^ + \Pt\ (3) 

We require that the transverse momentum contribution be less than the initial mass 
of the secondary vertex, to ensure that poorly measured vertices in uds events do 
not leak into the final sample by adding large Pt- 



5 Maximum Likelihood Fit 



Separation between the various lepton sources is accomplished using kinematic and 
vertexing information. Probabilities for each of the decay processes are assigned to 
every lepton, and calculated separately for electrons and muons. 



For electron candidates, eight discriminating variables are calculated based on char- 
acteristics of the event |^. These are track momentum (p), momentum transverse to 
the nearest jet (pt), same hemisphere vertex mass, same hemisphere vertex momen- 
tum, same hemisphere vertex significance, ratio of the track longitudinal distance 
from the IP along the vertex axis to the vertex distance from the IP (L/D, see fig. |l]), 
estimate of the underlying b quark boost, and the opposite hemisphere vertex mass 
(fig. 0). (Note: vertexing variables are not always all available for every event, but 
there is a requirement on the presence of a reconstructed vertex in at least one hemi- 
sphere.) These variables are fed into an Artificial Neural Network, configured with 
1 input layer, 2 hidden layers, and 4 output nodes (see figs. ||, ^ |10|, [TT|). The Neural 
Network weights are set by training on a Monte Carlo simulation of semileptonic 
decays of heavy quarks in Z° decays. The output of the Neural Network is checked 
with data. The probabilities are assigned by transforming the 3 NN output nodes 
onto a 2 dimensional space by x = A^A''^ + NN^c and y = NNi, + A^A"c, illustrated 
graphically by figs. |1^ and 0. This space is divided up into bins and every can- 
didate is assigned classification probabilites based on the number of events of each 
type in its corresponding bin. 



For the muon candidates, a multi-variate analysis is applied |^. Decay probabili- 
ties are calculated for every muon in the data by using a nearest neighbours technique 
in a 3-D Monte Carlo phase space. Three planes are defined, corresponding to three 
different ranges of the event mass (defined as the largest of the masses reconstructed 
in the two opposite hemispheres of the event, fig. §): > 2 GeV, 2 > M^tx > 0.55 
and 0.55 > M^tx > (or no vertex found) [|. These planes are parametrized by the 
quantities a/^ and In |-P|/2, to ensure a more uniform point distribution (also the 
scales for the two variables are roughly the same, see fig. ^. The weights for a muon 
in the data are then calculated with a nearest neighbours technique, by selecting 
all MC events within a fixed distance f\ from the data point in the corresponding 
plane and deriving the fractions of events of each type in this sample. In a second 
step and for those events only with a reconstructed vertex mass, these probabilities 
are re- weighted with fractions derived from the MC L/D distribution (see fig. ^, 
which account for the likelihood of an event coming from a certain source to have 
a value of L/D included in a specific interval. This information helps particularly 
to enhance the ability in separating b direct decays from b cascade decays (which 
have the "wrong" charge association). Correlations between all the different phys- 
ical quantities employed have been taken into account. A cross-check via a neural 
network approach (similar to the one used for the electrons) has been performed 
and it has given consistent results. 



A maximum likelihood analysis of all hadronic Z^ events containing leptons is 
used to determine Af,. The likelihood function contains the following probability 

^There is no vertex requirement for the muon analysis. 
^ (optimized with respect to the statistic available) 



term for each lepton in the data: 

P ip,pt, mass, L/D P„ z; A, A,) cx= {(1 + z'')il - A,P,) - 2Q{A, - P') 

\fb - he + fb,){l - 2xi){l - ^QCDi^))A, 
+/,(! - A'QCDiz))A, + ftkgAhkg] z} . (4) 

where z = cos 9 jet- The three signs governing the left-right forward-backward asym- 
metry — beam polarization Pe, lepton charge Q, and jet direction cos 6' — are 
incorporated automatically into the maximum likelihood probability function. The 
fractions (/&, /be, fbc, fc, hkg) are the lepton decay probabilities for the different de- 
cay modes, derived from the MC simulation as described before. A correction factor 
(1 — 2xi) is applied to all b-quark lepton sources to account for asymmetry dilution 
due to B^B^ mixing. For the electrons an average x = 0.1260 is calculated from 
Monte Carlo (to account for the analysis bias introduced by the vertex require- 
ment), and rescaled for h cascade events to account for different mixing probabilities 
(x = 0.1364 and x = 0.1376 are used for 6 — i> c — >■ e and 6 — >• c — >■ e events respec- 
tively). For the muon analysis, x is calculated event by event using the truth mixing 
information of MC events closest to the data point in a phase space parametrized 
by {p,pt, mass, L/D). The dependency on the event lifetime is thus automatically 
accounted for without any further need of rescaling. The background asymmetry 
Ahkg is derived for the electron analysis as a function of p and pt from tracks in the 
data not identified as leptons. For muons instead, it is calculated as a function of 
p and Pt (using the same procedure) from MC true background muons, divided in 
two samples: misidentified muons and light hadron muons (or misassociated tracks). 
A cos6 dependent QCD correction factor is applied to the theoretical asymmetry 
function to incorporate known QCD corrections to the cross section. The quantity 
Aq^j^{z) has been calculated at 0{as) for massive final state quarks by Stav and 
Olsen |2^ and is as large as 5% for the b quark at z = 0. For an unbiased sample 
of bb events with \z\ < 0.7, correcting for this effect increases the asymmetry by 
3% overall. However, the theoretical calculations assume perfect efficiency in the 
reconstruction of events with emission of gluons of any energy. The inefficiency of 
the detector and the presence of cuts and weights in the analysis cause biases in the 
event selection which favor qq events over qqg events, therefore the QCD correction 
to be applied is less than the theoretical one. The effects of this and other related 
biases have been studied with a MC simulation of the analysis chain and corrected 
for in the likelihood function as a function of 6, decreasing the theoretical QCD cor- 
rection by about 30% ||2^. 0{al) QCD effects are mainly due to two contributions: 
gluon splitting and second order hard gluon radiation. The gluon splitting correc- 
tion is calculated apart by re-fitting for Af, (in a MC as data study) after having 
excluded from the MC all the g —>■ bb, g ^ cc events. The difference in the central 
values, rescaled by the ratio of the current world measurements (OPAL) |]T| of the 
gluon splitting fractions to the JETSET input values, is assumed as correction. For 
the second order gluon radiation effects, recent theoretical calculations by Ravin- 
dran and van Neerven p8| have been implemented. These have been worked out 
for different values of the quark masses (pole masses or running) and they predict 



effects about four times as big as in previous calculations (~ 1% correction on Ai,) 



6 Results and Systematic Errors 



Tlie results obtained for tlie 1993-98 data are as follows, where the combined re- 
sult takes into account the systematic correlations between the muon and electron 
analyses. 



Muons 


A, 


Electrons 


A, 


Combined 


A, 



0.950 ± 0.038,t,t ± 0.026,^,i (5) 

0.876 ± 0.045,i,i ± 0.028,,,i 
0.922 ± 0.029,tat ± 0.024,y,i 



A list of systematic errors is shown in Table |T[ The background levels have been 
studied with the MC, but also with a data sample of pure pions from K'^ decays. 
The asymmetry of the background has been varied by ±40% of itself for the elec- 
tron analysis, and just rescaled by the ratio of the asymmetry in data and MC for 
charged non-leptonic tracks in the muon case. Uncertainty in the jet axis simulation 
can affect the asymmetry measurement by distorting the lepton pt spectrum and, 
to a lesser extent, the jet direction. The resulting systematic error has been studied 
by comparing the back-to-back direction of jets for data and MC in two jet events. 
The electron sample is more sensitive to such effects since both jet finding and elec- 
tron identification algorithm rely on the same calorimeter response. The precision 
of the 5^ and B^ lepton spectra is directly related to the uncertainty in the D** 
branching fraction reported by the CLEO collaboration |^|. The systematic error 



due to uncertainties in the D lepton spectrum has been estimated by constraining 



the ACCMM model [§ to the DELCO D -^ I data [g3|. The systematic error due 
to the QCD correction includes uncertainties in the 2nd order QCD calculations 
for hard gluon emission and gluon splitting, in the value of a^, and in the bias due 
to event selection criteria in the analysis. This analysis is independent of tracking 
efficiency, unless such efficiency depends on p, pt or is not symmetric in cos^. The 
extent of this p and pt dependence has been constrained by reweighting MC tracks 
by the ratio of the number of tracks in data and MC as a function of p and pt- The 
extracted value of Af is much less sensitive to potential differences in the relative 
efficiency for selecting leptons between the forward and backward hemispheres than 
are the values of Aj extracted from the unpolarized forward-backward asymmetry. 
The relative suppression factor is greater than 1/ A^^ ~ 50 for any value of |z| and 
therefore forward-backward asymmetry in the detector acceptance is not a signifi- 
cant source of measurement bias. A^ has been fixed in the maximum likelihood fit 
to its Standard Model value, and a systematic error has been calculated by varying 
this number by plus or minus twice the current statistical uncertainty on the world 



average of the Ac measurements. 



The value obtained for A^ from leptons can be combined with the other mea- 
surements performed at the SLC/SLD, respectively based on a momentum weighted 
track charge method, a vertex charge method and kaon decays. The resulting SLD 
average 



A = 0.914 ± 0.024 



obtained using the data collected in 1993-1998, is consistent with the SM prediction 
Ah = 0.935 and in agreement with recent preliminary results from LEP and SLD|l|. 



7 Conclusions 



In conclusion, we have measured the extent of parity violation in the coupling of 
Z^ bosons to b quarks by using identified charged leptons from semileptonic decays. 
The analysis presented in this paper is based on the entire sample of 550,000 Z'^ 
decays collected in 1993-98 at SLD and employs vertexing information to separate 
the different decay sources. The resulting 1993-98 measurement 

Ah = 0.922 ± 0.029 ± 0.024, 



represents an improvement relative to previous measurements 30 
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Tables and Figures 



Source 


Parameter variation 


5Ab{fi) 


5Ab{e) 


Monte Carlo weights 


fb, fc variation 


±.006 


±.006 


Track efficiency 


MC-data multiplicity match 


±.008 


±.001 


Jet axis simulation 


10 mrad smearing 


±.001 


±.010 


Background level 


±10%(/i), ±5%(e) 


±.003 


±.006 


Background asymmetry 


±40% 


^.003 


^.004 


Neural net training 


10 training runs 


±.000 


±.012 


BR(ZO -^ hi) 


Rb = .2173 ± .0007 


^.000 


^.000 


BR(ZO -^ cc) 


i?c = . 1674 ±.0038 


±.001 


±.001 


BR{b -^ I) 


(10.62 ±0.17)% 


^.004 


^.003 


BR(6 -^c^l) 


(8.07 ±0.25)% 


±.003 


±.003 


BR(6 -^c^l) 


(1.6 ±0.4)% 


±.005 


±.001 


BR(6 -^T^l) 


(0.452 ± 0.074)% 


±.002 


<.001 


BR(6 ^ J/iP -^ I) 


(0.07 ±0.02)% 


±.003 


±.002 


BR(c -^ I) 


(9.85 ± 0.32)% 


±.002 


±.002 


B lept. spect. - D** fr. 


(23 ± 10)%, B+,B^; (32 ± 10)%, Bs 


±.004 


±.003 


D lept. spect. 


ACCMMl iti^^c^^l) H 


±.004 


±.005 


S-tag 


eff. calibration 


±.014 


±.012 


L/D 


DT/MC ratio 


±.002 


±.000 


B^DD 


(7.2 ±2.0)% 


±.008 


±.000 


D^D^ 


15% uncertainty 


±0.000 


±.001 


Bs fraction in bb event 


.115 ± .050 


±.002 


±.005 


Afo fraction in bb event 


.072 ± .030 


±.002 


±.003 


6, c fragmentation 


Cb = .0045-.0075 


±.003 


±.002 


Aleph fragmentation 


ec = .045-.070 


^.003 


^.003 


Polarization 


<Pe>=. 729 ±.0038 


^.005 


^.006 


Second order QCD 


^QCD uncertainty 


±.005 


±.005 


gluon splitting 


Qbh, Qcc uncertainty 


±.002 


±.003 


B mixing x 


X= .1186 ±.0043 


±.015 


±.012 


Ac 


0.667 ± 0.040 


±.002 


±.005 


Total Systematic 




.026 


.028 



Table 1: Systematic errors for the maximum likelihood analysis (1993-98) 
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Figure 1: Topological parameters of a track: D is the distance of the secondary seed 
vertex from the interaction point along the line connecting them; T is the transverse 
distance of the track from the vertex axis calculated at the point of closest approach 
(POCA) and finally L is the distance from the IP of the projection of the POCA on 
the vertex axis. 
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Figure 2: Total momentum distribu- 
tion of identified muons in 1996-98 
data (dots) and Monte Carlo (his- 
togram). 



Figure 3: Transverse momentum dis- 
tribution of identified muons in 1996- 
98 data (dots) amd Monte Carlo (his- 
togram) . 
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97-8 Event Vertex Mass (zvtop3) 
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Figure 4: Event mass distribution for 
muons in 1996-8 data (dots) and Monte 
Carlo (histogram). The event mass is as- 
sumed to be the highest of the vertex 
masses found in the two hemispheres. 
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Figure 5: L/D distribution for muons in 
1996-98 data (dots) and Monte Carlo (his- 
togram) . 
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Figure 6: Distribution of 0.5 In |p|and ^/pr for niuons in the 1996-98 data 
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Figure 7: Vertex mass and L/D distributions for electrons in the data (dots) and 
Monte Carlo (histogram). 





Figure 8: Electron NN out- 
put, b ^ e node. 




Figure 9: Electron NN out- 
put, b -^ c ^ e node. 




Figure 10: Electron NN out- 
put, c — > e node. 



Figure 11: Electron NN out- 
put, background output. 
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Figure 12: Dalitz space plot of Neural Network output for data and Monte Carlo 
candidates. 
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Figure 13: Dalitz space plot for Monte Carlo candidates. Ideally, b direct decays 
should be clustered around (1,1), b cascade decays at (0,1), charm decays at (0,1) 
and background events at (0,0). 
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